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Abstract. 

Casimir friction between a polarizable particle and a semi-infinite space is a delicate 
physical phenomenon, as it concerns the interaction between a microscopic quantum 
particle and a semi-infinite reservoir. Not unexpectedly, results obtained in the past 
about the friction force obtained via different routes are sometimes, at least apparently, 
wildly different from each other. Recently, we considered the Casimir friction force for 
two dielectric semi-infinite plates moving parallel to each other [J. S. Hpye and I. 
Brevik, Eur. Phys. J. D 68, 61 (2014)], and managed to get essential agreement with 
results obtained by Pendry (1997), Volokitin and Persson (2007), and Barton (2011). 
Our method was based upon use of the Kubo formalism. In the present paper we focus 
on the interaction between a polarizable particle and a dielectric half-space again, and 
calculate the friction force using the same basic method as before. The new ingredient 
in the present analysis is that we take into account radiative damping, and derive 
the modifications thereof. Some comparisons are also made with works from others. 
Essential agreement with the results of Intravaia, Behunin, and Dalvit can also be 
achieved using the modification of the atomic polarizability by the metallic plate. 
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1. Introduction 

The Casimir friction problem - discussed extensively in recent years - is a delicate 
problem in fundamental physics. The typical geometrical setup is that of two parallel 
dielectric half-spaces moving longitudinally with respect to each other with constant 
velocity at close spacing. In the limit when the material density of one of the plates 
goes to zero, one obtains effectively the situation where one single particle moves parallel 
to a dielectric half-space. 

There are essentially two different avenues that one can follow to deal with the 
problem. One, probably the most conventional one, is to make use of field-theoretical 
methods. A second way is to make use of statistical mechanical methods from the 
beginning. This is the method that two of us have been using ourselves in several papers; 
it implies use of the Kubo response formula. In our opinion the statistical mechanical 
method is quite compact and useful in the Casimir context in general. In the present 
paper our main purpose is to extend our recent study for dielectric half-spaces DEI to 
the case where the influence from radiative damping is taken into account. However, the 
presence of the dielectric half-space is an additional source of energy dissipation that 
turns out to outweigh the radiative damping in free space by far at close separation. As 
mentioned below, this will be considered in two appendices. We limit ourselves to the 
case of a moving polarizable particle close to a half-space, and assume zero temperature, 

T = 0. Use of the effective polarizability for a harmonic oscillator will now be needed. 

The oscillator is damped by the radiation reaction due to the energy loss from radiated 
electromagnetic dipole radiation. Whereas the solution of this classical problem is well 
known, we mention that some time ago it was rederived by Hpye and Stcll within the 
framework of statistical mechanics for dielectric fluids [3j. It ought to be emphasized that 
in the present work the polarizable particle with its radiation damping will replace the 
influence from one of the metal slabs with its finite conductivity considered in Refs. H0. 

We work out the consequences from use of a complex polarizability. 

As results found in the literature are sometimes very different, there is a need for 
trying to make contact between the various approaches. We mention that we have earlier 
managed to get quite close agreement with results obtained earlier by Barton @H3, 

Pendry Pffl]], and Volokitin and Persson P3H35]. To our knowledge this has not been 
made before. As for the specific problem considered in the present paper, we have made 
contact with the recent calculation of Intravaia et al. m- The damping mechanism 
there is quite different, being due to the modification of the atomic polarizability by 
the conducting plate, and in fact in the two Appendices of this paper, we rederive their 
result based on this alternative mechanism. 

It may be of interest to give a few remarks about the historical development of this 
field. It is probably fair to say that the research on Casimir friction was started by the 
1978 paper of Teodorovich [18]. Our own research on Casimir friction was initiated early, 
by considering the simple but basic system consisting of a pair of harmonic oscillators 
that interacted weakly and moved slowly relative to each other. Both static [T9] and 
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time-dependent interactions [20] were considered. The first papers of Pendry from 
the 1990’s [8|[9j were influential, and established a kind of approach later followed by 
Volokitin and Persson, and others. A number of papers from the last 5 years are 
listed in Refs. [2T1 [37j . The research literature is quite diversified; some papers have 
focused on the complicated situation where the relative velocity is relativistic and even 
necessitating the inclusion of Cherenkov-like effects. We will, however, limit ourselves 
to the case where the velocity is nonrelativistic. Even this apparently simple case turns 
out to be intricate enough. 

The next section surveys the kind of statistical mechanical formalism that we have 
been developing in a series of papers. In section 3 we focus on the extension of the theory 
implying use of the radiation damping term (instead of finite conductivity for a plate) 
in the polarizability for the particle, pointing out how the conventional proportionality 
of v 3 in the zero temperature friction force becomes replaced by a proportionality of v 5 
(v is velocity). Some numerical evaluations and comparisons to other works are given 
in section 4. In Appendix A, we set up a field theoretic formulation, similar to that of 
Intravaia et al. mi. and use it to derive not only the radiation reaction result found here, 
but also the induced damping friction found in Ref. [17] (apart from a numerical factor). 
In Appendix B we rederive that same result using the statistical mechanical method 
explained in the body of this paper. Thus, there is essential agreement concerning the 
formulae for Casimir friction. 

2. Evaluation of friction 

For readability we shall first outline the main properties of the statistical-mechanical 
Kubo formalism that we are pursuing. The geometric setup is the conventional one: 
the lower plate 1 with surface situated at 2 = 0 is at rest, while the upper plate 2 with 
surface situated at z = d moves with nonrelativistic velocity v x = v in the x direction. 
We let the origin be located at the surface of the lower plate, and let r 0 denote the 
relative initial position between two particles in the plates. Thus r 0 = {xi,2/i,Zi — z 2 j, 
with x 2 = i /2 = 0. The friction force between a polarizable particle and a half-plane, as 
mentioned above, can be obtained when one half plane (chosen to be the upper one) is 
taken in the dilute limit. 

Let pi and p 2 be the particle densities in the two half-planes, and let kj_) 
be the Fourier transform in the xy-\Aaxie of the ^(r) ( r — r(i)) P ar t of the interaction 
-0(r)siS2 between a pair of oscillators. The kj^ is the Fourier variable, zq = Z\— z 2 , and 
Si and S 2 are the vibrational coordinates. The dissipated energy A E can be written in 
the form (cf. Eq. (23) in Ref. [l]) 



( 1 ) 


2l>d,Z2<0 
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As it can be shown that with electrostatic dipole interaction 

f2Tre- q \ z °\ \ 2 

kaMA), -k ± ) = (2 q 2 ) 2 (-j , (2) 

where q = |kj_| and oj v = k x v , k x being the component of the Fourier variable in the 
x direction, the z\ and Z 2 integrations of Eq. (JT]) can be performed. With this the 
dissipated energy becomes 

A E = Pl p 2 J e~ 2qd J{cu v ) dk ± . (3) 

Then consider the quantity J(u v ) in Eq. (JT|). The evaluation of it is more 
demanding, and involves use of the quantum mechanical response function for a pair of 
oscillators. Moreover, the quantity J(oj v ) contains information about the specific model 
that we use for the motion. In order to avoid other contributions to the change in energy, 
apart from dissipation, the motion should be performed in a closed loop. Without loss 
in generality we can assume that the motion takes place in the x direction. Initially, at 
t < —t — ar where r 1, a 1, we assume that the upper plate is at rest. In the 
period t G (— r — ar, —r), it moves to the left with a low constant velocity — v/a {y is 
constant and positive). Then, for t G (—r, +r), it moves to the right with finite velocity 
v , whereafter it moves to the left again with low velocity — v/a until it returns to its 
initial position at t = r + ar. In this way the closed loop motion of the upper plate is 
brought about. 

In the limit a —> 00 dissipation from the initial and final very low velocities can be 
neglected, and it turns out that we can write J(u> v ) as 

J(w„)=CJ(w.)+C + I(w + ), (4) 


where the coefficients C + and C _ are in general, at finite temperature, 


C± = — sinli 
n 


-/3hu± 


with u± = |aq ± u 2 \ and 


H = 


h 2 uiu)2aia2 


(5) 


4 sinh(i/3fta;i) 

Here a± and a 2 are the polarizabilities of a pair of particles, one in each half-plane, oj\ 
and ix >2 are the corresponding eigenfrequencies, and /3 = 1 /{ksT). 

As for the quantities I( 00 ) in Eq. (j4]), it turns out that they take the form 
,2 


c 0 . 


I(u) = 7rr—[5(uj — oo v ) + S(uj + cu„)] 


OJ 


( 6 ) 


for large r —» 00 . At T = 0 as assumed in this paper, C_ = 0 by which the J(lo v ) for a 
pair of particles will be 

J(uj v ) = C + ttt—[5(uj + - oj v ) + 5(o; + + u v )\, (7) 

cu + 


C + = -ftwiUJ2aia2. 


( 8 ) 
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Now consider higher densities of particles. Then, the polarizabilities are generally 
to be replaced by 

27r P« ->■ (9) 

where £ is the permittivity and p as before is the number density. As usual, we consider 
only local permittivity. The case of nonlocal permittivity (spatial dispersion) has also 
been studied and can sometimes be of significance; cf., for instance, Ref. [38] . 

In general the oscillators have a frequency distribution. Then we have the 
replacement 

OL a a Ia (m 2 a ) d(ml), m a = hw a , (a = 1, 2). (10) 


With integration elements d{m \) d(m |) = 4A 2 mim2 du\ du 2 this gives the following 
expression for J(co v ) (u; + = |ay|) 


\^v | 

J(oj v ) = 2TTr\uj v \h 3 j (jj\UJ2m\m20in{m 2 f)ai2(m\) duj\. 

0 

The frequency spectrum follows from the imaginary part of a or expression 
higher densities. It is obtained via Eqs. (27)-(29) of Ref. [36]. One can write 

a(K) = f(K 2 ), K = ihui. 

Then it can be shown that the function /(A' 2 ) satisfies the relation 

/(*W 


with 


(11) 

J) for 

( 12 ) 

(13) 

(14) 

5 K + 0 (15) 

where f = — iuj , the c o p is the plasma frequency, and v represents damping of plasma 
oscillations due to finite conductivity of the medium. The frequency spectrum that 
follows from the imaginary part of Eq. (j9j) is via Eqs. ([T2l) - (lT5li for small m (—> 0) given 
by 

hu 


K 2 + m 2 


ai(m 2 )m 2 = A/(— m 2 + iy), m = hu = — iK , 7 —> 0 + 

7 r 


With the Drude model for a metal the permittivity is given by 

u’ 2 

£ — 1 + 


cot 


= Dm 2 , D = 


pi^nhup ) 2 


(16) 


3. Radiative damping included 

After having presented this brief overview of our kind of statistical mechanical theory we 
now turn to the new element in the present paper, which is to include radiative damping 
for the polarizable particle moving above the half-space. (For a single particle finite 
conductivity and corresponding damping is obviously absent.) The damping results 
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from outgoing radiation from the particle. The resulting polarizability can be found, 
for instance, as Eq. (17) in the paper of Hpye and Stell on dielectric fluids |3], 

a 

a e = - 5 — 

l + ¥{f> a 

(note that Gaussian units are used). The ui in the reference is the —oj/c above where 
c is the velocity of light. (The minus sign is related to the different choice of sign 
by underlying Fourier transform.) It should be noted that this damping is for a 
polarizable particle in free space. However, it will be modified due to the presence 
of a dielectric medium where energy can dissipate. This extension is considered in two 
appendices where different methods are applied. For close separation the influence from 
the dielectric medium turns out to dominate by far. In Ref. [3] expression m was given 
a statistical mechanical derivation, but it is consistent with the one obtained via the 
radiated power from an oscillating dipole. Here a is the polarizability of the undamped 
oscillator 

a = —y -; 

(Uq — UJ 

For small frequencies (lo -C co 0 ) needed here (v c) the a can be replaced by a 0 by 

which Eq. (1171) becomes 




ot e —y oil 


a 


i+thtr 


— a o 


«o 




ol o 

(he ) 3 



/(-m 2 +i7)(19) 


With Eqs. (TT2li and ifldli this gives the frequency distribution (7 —> 0 +) 

o; n (m 2 )m 2 = --$/(—)m 2 + 17 ) = Bm\, B = (20) 

7 r o7T[nc)° 

Expressions (OR)) and (12U|) for the frequency distributions are to be inserted in 
integral (fTTlh Inserting Eqs. (TT6l) and (120|) in Eq. CD we obtain ( uj 2 = |a 7 ,| — u>i) 


\(Bv | 

J(u v ) — 2nT\ui v \h 5 BD I ojloj 2 doJi = 2 tu^H pii , H PII = -— h 5 BD. ( 21) 

J 20 

0 

For clarity, an extra subscript II is given to quantities differing from corresponding 
quantities calculated in Ref. p]. The most notable difference from Eq. (52) in |lj is 
that Eq. (1211) contains the factor cu® instead of Cod (the factor r present in H P in [I] 
is a misprint). Note that in Ref. [Tj two equal metal half-planes were considered by 
which Eq. (1T61) was used for both while in the present case Eq. (|20l) is used for one 
of them. By further integration of cu® = (k x v) e (v = v x ) over orientations one has 
IK df = k* ± J 0 27T cos & (pd(j) = 27rg 6 (5/16) ( q = k±) such that k x can be replaced by 
5g 6 /16, i.e. we can write 

JM = 2rv % H PI Xq\ ( 22 ) 

lo 

which replaces Eq. (53) in [I]. The most characteristic difference is that v 6 q 6 replaces 
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We can now calculate the dissipated energy per unit area from Eq. 
AEpn = 2rv 6 HpuG pij, 
where Hpn is given by Eq. (j2Tj) above, and 


to be 


OO 

Gpii = P 1 P 2 j yjr<? 6 e~ 2qd 2vr q dq = p x p 2 
o 

The friction force per unit area becomes accordingly 

AEpu 


5 • 315tt 
2 7 d 8 ' 


Fpn — 


— —HpnGpuv — — pip 2 


315tt 2 h 5 


■ BDv 5 


(23) 


(24) 


(25) 


2 TV rir ‘ 2 9 d 8 

The results above were for two half-spaces. Consistent with the result (125]) the 
force between a particle moving at a distance z 0 away from a half-space will be given by 
(extra subscript II omitted from now on) 


^=-4 

^0 


where the coefficient A follows from Eq. (1251) as 


3157r 2 h 5 5 105hz/«o 5 

A = 8 p 2 ——— BDv 5 = — - p-Z-v 5 


2 9 


2 5 7 TUpC 3 


(26) 


(27) 


(We have here made use of the fact that the force for two half-spaces results from 
integrating the single-particle-half-space force over all z 0 , from d to infinity.) Recall that 
B and D are given by Eqs. (1151) and (12U[) respectively. Note also that these expressions 
are given in Gaussian units, implying that the polarizability a 0 has the dimension cm 3 . 

We emphasize the difference between the present result (1251) . oc v 5 , and the earlier 
result oc v 3 obtained in [T], Eq. (56), for two metal half-planes with use of the same 
dielectric function (1T5]) . The reason is the polarizability (fT7|) for the low density medium. 
With its (n 3 -behavior in the denominator this differs in a qualitative way from expression 
(1I5|) that gives a corresponding cn-behavior in expression ([9]). This has the consequence 
that the frequency distribution has the Bm\ behavior of Eq. (1201) while the latter has the 
Dm 2 one of Eq. (116]) (m = hut). This difference in behavior reflects the rapid vanishing 
of the radiation damping as the frequency decreases. 


4. Remarks on numerical magnitudes 

Various approaches to the Casimir friction problem sometimes look very different from 
each other. This may not be unexpected. Due to this it may be difficult to check 
whether disagreements between results obtained are real or merely reflect differences in 
formalism. Little effort has to our knowledge been made by others to clear these things 
up. As mentioned, our method is based upon the Kubo formalism [40] . 

To begin with, let us present some encouraging results for the case of two half¬ 
spaces. Volokitin and Persson, in their 2007 review article [13], found the T = 0 friction 
force here called Eyp, to be Fyp = F/2, where F is our result [1,36]. The expression 
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Fyp was based upon a calculation of Pendry [ 8 ], who found Fp e ndry = F VP /6 (this slight 
mismatch is seemingly due to a trivial calculational error in Pendry’s paper). The same 
situation was considered also by Barton [7], who found Fhai-ton = 12Fp en dr y (when a 
factor ((5) = 1.037 is disregarded). So, we find actually agreement at T — 0, 

F = Fhai-ton (two half spaces). (28) 

Then turn to the case which is more relevant to the present case, namely one 
particle moving outside a half-space. It would here be of interest to make contact 
with the numerical results obtained recently by Intravaia et al. ini. who considered a 
situation of the same kind. Now, it is known that the largest polarizabilities for the 
elements are found for alkali atoms. We choose rubidium as an example (the same 
element as considered in Ref. 031)- The static ground-state polarizability op ( 5 2 S'i/ 2 ) is 
in this case given in Ref. [3TJ as 

a 0 = h x 0.0794 Hz/(V/cm) 2 , (29) 

with h = 6.626 x 1CU 34 Js being Planck’s constant. Thus, in ordinary SI units, 

a 0 = 5.26 x lO ” 39 Fm 2 , (30) 

where F=C/V means farad (the dimension given in Ref. [TTJ is incorrect). Using the 
conversion formula between Gaussian polarizability ao —> «og and SI polarizability ao 

«og = t -atb (31) 

47T£o 

with the permittivity of vacuum e 0 = 8.85 • 10 -12 F/m, we get the Gaussian ground-state 
polarizability 

a 0G = 47.3 x 1 (T 24 cm 3 = 47.3 A 3 , (32) 

in terms of the convenient unit 1 A = 1 CT 8 cm. (Most other elements have a static 
polarizability of order 1 A or less.) 

As for the plasma frequency and the relaxation parameter we choose for definiteness 
the values for gold, hio p = 9.0 eV and hu = 35 meV or u> p = 1.36 x 10 16 rad/s and 
v = 5.32 x 10 13 rad/s (cf., for instance, Ref. [47]). Then, choosing v = 340 m/s we 
obtain from Eq. (17711 

A = 1.19 x 1(T 122 Nm 9 = -1.19 x 10~" dyne cm 9 . (33) 

Further choosing zq = 10 nm we obtain from Eq. (|26l) the force on a single particle to 
be 


F = -1.19 x 10" 5 ° N = -1.19 x 10 “ 45 dynes. (34) 

This is far beyond detectability in practice, and it is far lower than the force found 
in Ref. [IT]. However, the latter result was found for a silicon half-space instead of 
gold. For silicon the resistivity p = 6.4 • 10 2 Urn was used, and with the other numbers 
unchanged the friction force F zz —1.3 • 10~ 2 °N was found. Likewise for the same 
situation with use of the relation p = u/(£ 0 ojp) our Eq. (1261) together with (1771) give 
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the tiny friction force F = —2.3 • 10 -40 N. Although larger than result (EHj) it is still 
far from the result of Ref. m- However, as shown in the Appendix A, our result can 
also be derived using the formalism of Ref. m, so we conclude that the difference 
resides in the different atomic dissipation mechanism assumed, and indeed, by using the 
modification of the polarization induced by the metallic plate, we recover (essentially) 
the result of that reference as well. This is verified by a calculation in Appendix B which 
again reproduces (essentially) the result of Ref. [T7] using the formalism given above. 

Thus the controversy appears resolved; in practice, the mechanism based on the induced 
particle dissipation by the interaction with the plate is by far larger than that caused 
by radiation reaction, although below the level of detectability with present techniques. 

We may conclude our work as follows: Our intention has been to evaluate the 
friction force at T = 0 for a polarizable particle moving parallel to a metal half-space 
at close separation, using the statistical mechanical Kubo formalism. First, friction 
connected with the radiation reaction in free space is found. Then, the influence 
from dissipation in the metal half-space is included in the appendices where also an 
independent method of derivation is used. At close separation the latter dissipation is 
found to outweigh the former by far. Moreover, the results obtained by our method are 
compared with results obtained by others, based upon independent and very different 
methods. Essential agreement is obtained. 


Appendix A: Field theory (source theory) approach 


In this appendix we wish to make contact with the formalism employed by Intravaia et 
al. HT] and rederive the above result for Casimir friction between a polarizable atom 
and a metal plate. Because the approach seems quite different, and the derivation in 
Ref. [17] is quite terse, we describe it in detail. (Here we set h = c = 1 and use, except 
where otherwise noted, rationalized Heaviside-Lorentz units.) 

We start from the action term describing the interaction between two polarization 
sources, 


W = (dx)(dx')P(x) • T(x,x') • P(z'), (A.l) 

where the integration is over all space-time coordinates, and T is the Green’s dyadic 
for the system. The atom offers a realization of the source P; if the atom moves with 
velocity v, as before, parallel to the surface, the change in the polarization in a small 
time 5t is 5P = —5tv ■ VP. The resulting variation in the action allows us to read off 
the force due to the quantum friction in terms of the power (T is the “infinite” time 
interval the system exists) 


P = F • v 


(HR 

T8~f 


(A.2) 


where 


F = 


/ 


d (t — t , )(dr_ L )(dr( L )d 2 ; dPP(rj_, z, t ) • (V)r(rj_, z; r' ± , z'\ t — t') ■ P(r'j_, z\ t').( A.3) 



Casimir friction between polarizable particle and half-space with radiation and image damping at zero temper 


Now we do a 2 + 1 dimensional break-np of the Green’s dyadic for the surface, 
T(Y L ,z]v' L ,z'-t-t') = - 1 - [ due-“M [ l£^j* x < rx -^ ) g(z,z'ik ± ,u), (A.4) 


2tt 


(2tt) 


and we replace the polarization by its realization by a moving dipole, 

P(r ± , z, t) = d(t)S(r ± - r ± (t))d(z - z 0 ), (A.5) 

where the atom is a distance z 0 above the plate. Here r_|_(t) is the atom’s trajectory, 
which we approximate by motion with constant velocity, r ±(t) = v±t, in which case we 
obtain the simple formula for the frictional force 


F = 


den 

27T 


(dkx). 


(2n) 


ik 


l 


dre -iG- k X-vx)r tr[d(i / )d(t)g(2o , 2o) ], (A _ 6) 


with r = t — t'. We interpret the product of dipole moments here as a correlation 
function, 


C(r) = -(d(f)d(f') + d(Od(f)), 

and we see the appearance of its Fourier transform, in the notation of Ref. m 

i 


S(w) = 


2 vr 


dr e iajr C(r). 


Thus our final formula nearly coincides with that given in Ref. 

(dkj 


F = 


den 


(2tt) 


-ikj_ tr[S(en - k ± • v ± ) ■ g(n 0 , z 0 ; k ± , en)]. 


(A.7) 


(A.8) 


(A.9) 


The power spectrum S is related to a generalized susceptibility through the 
fluctuation-dissipation theorem (FDT). We follow Landau and Lifshitz [43.J and write 
at zero temperature 
1 


S(en) = sgn(en)—Aa(eu), andsimilarly g(en) —> sgn(en)g(en). 
2ir 


(A.10) 

The formula 1IA.9|) differs in two crucial aspects from that of Ref. HZ|: The integral 
ranges over both positive and negative frequencies, and the sign of the argument of the 
power spectrum S is reversed. These modifications are necessary to recover the correct 
results. In fact, if we use (IA.9P to rederive the normal Casimir-Polder force for v = 0, 
by replacing ik± —)■ d z , we see immediately that the expected result is found, 


F z = —d z Ucp, U C P = 


dCtra G (i|C|)r(z 0 ,z 0 ;i|CI), 


(A.ll) 


where we recognize that the relation between the Heaviside-Lorentz (used elsewhere in 
this Appendix) and the Gaussian values of the polarizability is ckhl = 4770:^- All that 
is needed to establish this identity is the Kramer-Kronig dispersion relation 


«(*|C|) = - 

vr 


do; 


cuAa(a;) 


(A. 12) 


w 2 + C 2 

Then the frictional force can be shown to arise only from the low frequency region 
because the k^ integrand must be even: 

(dkj.) /‘ k v day 


F = 


(2tt) 


7T 


-k_L tr[AcK(a; — k • v)Ag(a;, fcj_)]. 


(A.13) 
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This result differs from the formula given in Ref. fa by a factor of and the sign of 
the argument of ct. 

We now assume that the velocity of the atom is very small, 

uj < k_L • v_l <C ujp, (A. 14) 

so that only the transverse magnetic (TM) part of the Green’s dyadic contributes, and 
for an isotropic atom we encounter only 


where 


trSg(z 0 j Zo; &,<*>) 
£(u) - 1 




£ — 1 

. -( 

£ + 1 


-2k z 0 


A 


2uv 


LOi 


(A.15) 


(A.16) 


£(u) + 1 

This says that the dominant effects of Casimir friction occur in the nonretarded 
evanescent regime. 

All of this seems completely unambiguous, and noncontroversial. It is the nature 
of the susceptibility of the atom that that seems to cause some consternation. The 
assumption being made in the body of this paper is that the dissipation mechanism in 
the atom is that of radiation reaction, so as noted in the text 


rv. / \ ^23 

= —Q'qOJ . 
o 


(A.17) 


Then when we substitute (lA.lOll . (1A.17I) . (1A.15I) . and (1A. 16|) into (1A.9D . and carry out 
the elementary integrations, we obtain (a factor of 47 t comes from using unrationalized 
Gaussian units) 


F = 


105 a 2 0 vv 5 
327T 0J% Zq 


(A.18) 


This coincides with the result (J2TJ) found in the text, and indeed the integrals involved 
are the same as encountered there. 

Quite different assumptions are made for the atomic susceptiblity in Ref. ia, 
leading to a very different dependence on atomic velocity and distance, F oc v 3 /d 10 . 
In fact, it is easy to reproduce that result. The point is that the image of the charge 
distribution in the atom lags behind that of the moving atom due to resistance in the 
metal. This is captured by rewriting (lA.ip symbolically as 


W = ^ Tr E • a 0 • T 0 • a 0 • E, 


(A.19) 


where a 0 is the polarizability of the atom in empty space, and T 0 is the Green’s dyadic 
of the metal plate without the atom present. Both are modified in the presence of each 
other, and then the effective source term is 

W e s = i Tr nT. yT = (EE). (A.20) 

Therefore, the effective polarizability of the atom in the presence of the plate is 


a e ff — ~2 a ° ' ' ao ’ 


(A.21) 
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or in the evanescent region, using (1A.16I) . 


SO! eff = 


On 


uv 1 


diag(l, 1, 2), 


(A.22) 


vr (2^ 0 ) 3 a;2 2 

where we have assumed that the atom is isotropic and possesses a real polarizability. 
The anisotropy arises from that of the TM Green’s dyadic. 

Now we insert this into the general formula (IA.9[) with (1A.10D for the power 
spectrum, together with the Green’s dyadic properties on the principal diagonal 


%g(z 0 ,z 0 ;k,uj) « ye 2fc2 °ydiag(l, 1, 2), 

2 £ + 1 2 


again with ( 1A.16I) . which immediately yields 


F = 


135q;qU 3 


(A.23) 


(A.24) 


4:7T 3 a 2 (2zo) 10 ’ 

which is 3/8 times the result of Ref. [IT]: This factor comes from the trace over the 
product of Green’s dyadics. Here we use the connection between the conductivity of the 
metal and the dissipation frequency in rationalized units, 

,2 


C 0, 


a = 


JL 

v 


(A.25) 


Thus the controversy involving the friction between an atom and a metal surface 
seems resolved. The different results depend on the damping mechanism assumed. In 
practice, it would in fact appear that the last mechanism, found by Ref. HU, and 
confirmed, again up to a factor, by Ref. is dominant for low velocities, although it 
still seems far beyond experimental reach. 


Appendix B: Statistical mechanics approach 


We also want to make contact with the result of Ref. [IT] using the methods of Secs. [2] 
and |3] With radiation reaction the method of Appendix A gave result (1A.18I) which 
fully agrees with results ([26]) and f[2T|) . 

The radiation reaction is an interaction on an oscillating dipole back on itself. 
In Ref. [3] it was found by statistical mechanical study of the refractive index of 
polarizable fluids that the radiating dipole interaction gave a self-interaction that adds 
to the oscillator potential by which the given polarizability is replaced by an effective 
polarizability. 

In the present situation there is a dielectric half-plane near the polarizable particle. 
This clearly absorbs radiation and thus adds to the damping of the oscillator. To obtain 
this damping we need the field that the half-plane reflects back to the oscillator. For 
small separations it turns out that this field by far outweighs the usual radiation reaction. 
Thus the latter can be neglected, and it is sufficient to consider the electrostatic limit 
of the electric field on a particle near near a dielectric half-plane. As noted above 
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Eq. (IA. 19[) the solution of this problem (outside the half-plane) is equivalent to adding 
a mirror oscillator in the half-plane. The electric held from a dipole moment s' is 


E = —-(3(s' •?)?-§'), 


(A.26) 


where the carets denote unit vectors. Let s' be the mirror dipole of the dipole moment 
s of the polarizable particle. For the components parallel and transverse to r we thus 
have (mirror picture has opposite charges) 

e — l 

Sj_ 


s II — crs ii 


with a = 


£ + 1 


With this the electric held is equivalent to a self-interaction 


0 =-J E.ds = ^ + a ; + 2 ^), 

with s z = S||, s\ = s 2 x + s 2 , and where r = 2zq. 

The oscillator with given polarizability «o oscillates in a potential 
1 


(A.27) 


(A.28) 


0o — 


2a, 


/ 2 i 2 , 2\ 

fs x + s + s z ). 


0 


The resulting potential is thus 0 O + 0 with effective polarizability (i = x, y, z) 
11 a 

®-ei 


ao 


(2z 0 y 

o <y 




with 


Tlx Tly 1 ■ 


n z = 2. 


&ei O'q a 0 - 


jUi -\ -. 


(A.29) 

(A.30) 

,(2^)3- ■ • ( A - 31 ) 

By the statistical mechanical study only the situation with scalar polarizability was 
considered for simplicity. In the present case where the induced part is anisotropic, 
one will expect an effective scalar polarizability to be equivalent to some average of it. 
The approach of Appendix A with anisotropy included tells how this average should be 
performed. According to the Green’s dyadic properties of Eq. (jA.23j) double weight is 
put on the ^-component compared to the x- and ^-components, i.e. the ^-component 
contributes most to the friction as it does to the induced polarizability. Thus with 
weighted average (n) = (n x + n y + 2n z )/A = 3/2 the average of a e i is 

o cr 


Ol e ~ OlQ 


Gin 


(n) + 


°(2.~ 0 )3 

The frequency spectrum is obtained like Eqs. (TT9l) and 

(n) = Dim i, D 1 = 


as (mi = hio) 


an(mi)mj = - 

71 


Qtn 


( 2 zo)' 


r Act 


3«q hu 


(A.32) 


(A.33) 


87T^o(^p ) 2 ’ 

where Eq. (1A.16|) is used. This frequency distribution is linear in mi like the one given 
by Eq. (ITfil) for the half-plane. By evaluations similar to Eqs. (T2Tji — (1271) one finds the 
friction force to be given by Eq. (56) of Ref. [Tj with one factor pD replaced by pi-Di 
where p\ —> 0 is the low density of polarizable particles forming the other half-plane 
utilized in Secs. [2] and [3] 

15?r piDipD{hvf 


F d = 


64 d 6 ‘ 


(A.34) 
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with D given by Eq. (fIBl) . Consistent with this and similar to Eqs. (126[) - (I27D the friction 
force on a single particle becomes 


F — — 


A 


with A = 6 • 


= S&fa) 


(A.35) 


2 8 ir(huip) 4: ZQ 

(By the previous derivation the Fn followed from F with D\ = D constant. Integration 
with varying Di would only change (1A.34D by a factor 6/9 while (1A.35D would remain 
unchanged.) This coincides with the result (1A.24I) with (I A.2 51) and conversion formula 
(l3Tf inserted (with units where £ 0 = 1 as used in Appendix A). Thus the result (up to 
the trace factor mentioned) of Ref. jTTj has been established in two rather independent 
ways where the arguments and reasonings are somewhat different. 
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